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X-ray Shadowing Experiments Toward Infrared Dark Clouds 

L. D. Anderson 1 , S. L. Snowden 2 k, T. M. Bania 3 * 

ABSTRACT 

We searched for X-ray shadowing toward two infrared dark clouds (IRDCs) 
using the MOS detectors on XMM-Newton to learn about the Galactic distribu- 
tion of X-ray emitting plasma. IRDCs make ideal X-ray shadowing targets of 
| keV photons due to their high column densities, relatively large angular sizes, 
and known kinematic distances. Here we focus on two clouds near 30° Galactic 
longitude at distances of 2 and 5 kpc from the Sun. We derive the foreground 
and background column densities of molecular and atomic gas in the direction of 
the clouds. We find that the | keV emission must be distributed throughout the 
Galactic disk. It is therefore linked to the structure of the cooler material of the 
ISM, and to the birth of stars. 

The diffuse X-ray background of the Milky Way has been studied for over 30 years and 
yet relatively little is known about the distribution of its emission in the Galactic plane. Most 
diffuse Galactic emission arises from ~ 1— 3x 10 6 K plasma, but the origin and distribution of 
this plasma are still debated. Two factors are responsible for the difficulty in determining the 
distribution of this plasma: confusion and absorption. Confusion arises because one cannot 
determine the origin of an individual X-ray. Even though most X-rays contributing to the 
diffuse background of the Milky Way are thermal in origin and the spectra of this plasma 
show T strong line emission, the available X-ray detector energy resolution is insufficient to 
provide velocity (and therefore distance) information. Thus it is very difficult to determine 
the true distribution of X-ray emitting plasma along a given line of sight. 

Shadowing experiments are the only way to determine the distribution of hot plasma 
in the plane of the Milky Way. To date, they have been used primarily in the | keV band 
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to study the distribution of 0.1 keV plasma within the Local Hot Bubble and in the lower 
halo (e.g., Snowden et al. 2000). A few shadowing experiments at higher energy (f keV 
and 1.5 keV) using ROSAT data (e.g., Park et al. 1997; Almy et al. 2000) have shown 
evidence for the existence of an extensive distribution of hot plasma well within the solar 
circle. This emission has been linked to the Milky Way X-ray bulge which has a scale height 
of ~ 1.8 kpc and a radial extent of ~ 5 kpc (Snowden et al. 1997), but there are likely 
additional components distributed throughout the Galactic disk. 

Most studies of the X-ray background have focused on the Local Hot Bubble (LHB). 
The LHB is a region of about 100 pc in diameter, surrounding the Sun filled with a hot, 
rarefied 10 6 K plasma (e.g., Snowden et al. 1990). Beginning with Snowden, McCammon, k 
Verter (1993), many groups have effectively mapped the distribution of this gas (e.g. Smith 
et al. 2005; Galeazzi et al. 2007; Henley, Shelton, k Kuntz 2007; Smith et al. 2007; Henley k 
Shelton 2008). These studies have shown the utility of X-ray shadowing in mapping a large 
region of the sky at higher Galactic latitudes. Here we utilize similar techniques but extend 
the analysis to larger distances in the Galactic plane using dense molecular clouds. 

Infrared Dark Clouds (IRDCs) are very dense molecular clouds that appear dark at 
infrared wavelengths. They are seen in absorption against the bright mid-infrared Galactic 
background. IRDCs were first identified as a significant population in the Midcourse Space 
Experiment ( MSX ) mid-infrared data (Egan et al. 1998), and are a ubiquitous feature of the 
higher resolution Galactic Legacy Infrared Survey Extraordinaire (GLIMPSE: Benjamin et 
al. 2003). IRDCs have very high column densities of 10 22 cm -2 (Simon et al. 2006, hereafter 
S06), and possibly as high as 10 23 — 10 25 cm 2 (Egan et al. 1998; Carey et al. 1998). They 
exist at a wide range of distances and are most easily seen in the plane of the Galaxy (S06). 

S06 found distances to 313 IRDCs using the Galactic Ring Survey of 13 CO emission 
(GRS; Jackson et al. 2006) to establish a morphological match between the 13 CO gas and 
the mid-infrared extinction. Using the rotation curve of Clemens (1985), S06 converted 
the velocity of the associated 13 CO gas into a distance to the IRDC. In the inner Galaxy 
(the focus of the present study), every velocity has two possible distance solutions. IRDCs, 
however, are seen in absorption against the Galactic plane and therefore the near distance 
can be assumed. 

Infrared Dark Clouds (IRDCs) make ideal X-ray shadowing candidates. Their high col- 
umn densities make them effective absorbers of background X-ray photons; X-ray emission 
from the direction of the IRDC must originate in the foreground. Thus, using IRDC ab- 
sorption, one can separate foreground and background X-ray components. The distribution 
of IRDCs also spans a wide range of distances (see S06. Jackson et al. 2008). With enough 
X-ray absorption measurements of IRDCs along different lines of sight, an accurate map of 



the hot ISM may be created. 


1. IRDC Sample 

The present study includes two clouds from the catalog of S06 that have high column 
densities, large angular sizes, and are closely grouped on the sky. Our goal is to select clouds 
that would cause complete absorption of X-ray photons originating behind the cloud. The 
models given in Snowden et al. (1994) predict that a total column density of 1 x 10“ 2 cm 2 will 
absorb ~ 98% of incident photons at | keV. Column densities of 10 22 are frequenctly found 
for IRDCs (see S06). To ensure sufficient counts in both the on- and off-cloud directions, 
the target IRDCs should cover approximately half the XMM-Newton detector. Our ideal 
situation is one in which there are multiple large clouds closely grouped on the sky but at 
various distances. The close grouping allows one to disentangle the emission components 
along a given line of sight. 

The properties of the two IRDCs observed by XMM-Newton are summarized in Table 
1. Listed are the IRDC name, the Galactic longitude and latitude, the angular size, and 
the molecular column density. All parameters in Table 1 are reproduced from S06. The two 
target clouds have high column densities, large angular sizes, and are near l = 30°. 

The mid- infrared emission of the two IRDC fields is shown in Figure 1. The background 
image in this figure is from the GLIMPSE 8/im survey (Benjamin et al. 2003) and the con- 
tours are 13 CO integrated intensity from the GRS. We have created the integrated intensity 
map using the LSR velocity and line width from the S06 catalog. Tick-marks on the contours 
indicate the direction of decreasing 13 CO integrated intensity. The large circle in this figure 
is approximately the XMM-Newton field of view (~ 30') . Immediately evident in this figure 
is the excellent morphological match between the GRS 13 CO contours and the strong mid- 
infrared extinction. The molecular clouds associated with the IRDCs take up a significant 
fraction of the XMM-Newton field of view. IRDC G28.37+0.07 is located at the center of a 
large, dense molecular cloud which is roughly ellipsoidal in projection. The tw r o clouds are 
at very different distances, G28.37+0.07 lies ~ 5 kpc from the Sun and G36.67— 0.11 lies 
~ 2 kpc from the Sun (S06). 


2. Observations and Data Reduction 

Our clouds were observed by XMM-Newton on 2006 April 11 for 25 ks each (IRDC 
G28.37+0.07 with XMM-Newton ObsID 0302970301 and IRDC G36.67-0.ll with ObsID 
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0302970201). 

One of us (SLS) has developed a powerful suite of PERL and FORTRAN scripts to 
analyze X MM- Newton extended source data called XMM-ESAS (Snowden et al. 2008). This 
software is based on the background modeling described in Snowden, Collier & Kuntz (2004) 
and Kuntz h Snowden (2008). The software currently only operates on the MOS detector 
data and subsequent to this analysis has been incorporated into the XMM-Newton mission 
Standard Analysis Software (SAS) package. 

To remove periods of high soft proton contamination, we must first temporally filter 
and clean the events file. The signature of soft proton contamination is a fluctuating light 
curve. The XMM-ESAS software creates a count rate histogram and fits a Gaussian to 
the peak of the histogram. In observations with minimal contamination, this histogram 
would appear Gaussian: a non-Gaussian shape to the distribution is also indivative of soft- 
proton contamination. The software defines count-rate intervals within 1.5er of the median 
count rate as acceptable, and filters the events file to remove other periods. While effective 
at removing intermittant contamination, this process does not remove consistent low-level 
contamination that is also likely present. The light-curve for G28.37+0.07 is quite steady 
and almost all the data are useable. The light-curve of G36.67— 0.11, however, indicates a 
high level of soft proton contamination and only ~ 20% of the observation time is usable. 


3. Estimation of Column Density 

An estimate of the total column density to each IRDC can drastically improve the 
reliability and accuracy of the X-ray spectral fitting. Above 0.5 keV, the most significant 
source of X-ray absorption is metals. With normal abundances, the total opacity should be 
proportional to N( H) + 2A T (H2). Thus, we estimate the molecular column density using the 
13 CO emission from the GRS, and estimate the atomic column density using the 21 cm H I 
VLA Galactic Plane Survey (VGPS: Stil et al. 2006). 

Our spectral fitting procedure treats the on- and off-cloud spectra separately. This is 
the only way one can separate components of the X-ray background that originate in front of 
the cloud from those that originate behind it. One must therefore compute on- and off-cloud 
column densities separately. We do this by first defining the on- and off-cloud regions. 

One of us (LDA) has developed an astronomical software program in IDL. 1 This soft- 
ware allows one to define regions of interest in an image using a threshold selection tool. 


1 Download at http : / / www. bn .edu/ iar / kang/ . 



The user can vary the threshold upwards or downwards to define smaller or larger regions. 
The output of this thresholding selection, the x and y positions of the data points on the 
threshold boundary in sky coordinates, can then be converted to detector coordinates within 
the software, as is appropriate for analysis in the XMM-Newton Science Analysis Software 
(SAS 2 ). This conversion employs the SAS program edet2sky repeatedly for each coordinate 
pair to convert from a region in sky coordinates to a region in detector coordinates. 

Using our IDL software, we define the on-cloud region as the sky positions possessing 
both high mid-infrared extinction and high 13 CO column density. Our goal is to create a re- 
gion large enough to have sufficient counts in the on-cloud region for good statistics. Ideally, 
the on- and off-cloud regions would each cover half the XMM-Newton field of view. The 
off-cloud region encompasses all locations in the XMM-Newton field of view not including 
the on-cloud region. The on-cloud region roughly corresponds to the largest contour in 
Figure 1. 


3.1. Molecular Column Density 

We must separate the total molecular column density into a component forground to 
the IRDC and a component beyond the IRDC. We therefore must estimate the distances to 
the molecular gas clumps that contribute to the column density in the XMM-Newton field 
of view. In the first Galactic quadrant (where the target IRDCs lie), each radial velocity 
has two possible distance solutions, a “near” and a “far” distance. This problem is known 
as the kinematic distance ambiguity (KDA). Gas along the same line of sight as an IRDC 
could lie either in front of or behind the cloud. This problem is illustrated in Figure 2 for 
IRDC G28.37+0.07. Figure 2 shows the loci of distances from the Sun corresponding to a 
particular velocity, assuming the Galactic rotation curve of Clemens (1985). The two- fold 
distance degeneracy is clearly shown. Because IRDCs are seen in absorption, we assume 
they lie at the near distance. The cloud velocity of 78.6 kms -1 places G28.37+0.07 at a 
distance of 4.8 kpc from the Sun. Any gas with a velocity greater than 78.6 kms -1 must 
lie beyond the cloud. Molecular gas at velocities less than 78.6 kms -1 , however, could lie 
at either the near or the far distance. 

Hi self-absorption (Hi SA) can remove the distance ambiguity of the molecular gas 
clumps in the XMM-Newton field of view. This technique relies on cold foreground H I ab- 
sorbing the emission of warm background Hi at the same velocity. Warm H I is ubiquitous 
in our Galaxy and emits at all allowed velocities. Litzt, Burton k Bania (1981) hypothe- 


2 http://xmm.esac.esa,int/sas/ 
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sized that molecular clouds must contain H I, a result which has been confirmed in numerous 
subsequent observations (e.g., Kuchar & Bania 1993; Williams & Maddalena 1996). This 
small population of cold neutral H I atoms inside molecular clouds is maintained by inter- 
actions with cosmic rays. The H I inside molecular clouds at the near distance will produce 
an absorption signal because there is ample warm H I emitting at the same velocity at the 
far distance. Any cloud at the far distance should not show Hi SA because there is no 
background H I emitting at the same velocity. H I SA is therefore a general technique for 
finding distances to molecular clouds. This technique was shown by Jackson et al. (2002) to 
be effective at resolving the KDA for a dense molecular cloud, by Busfield et al. (2006) for 
resolving the KDA for young stellar objects, and by Anderson & Bania (2009) for resolving 
the KDA for the molecular gas associated with Hu regions, and by Roman-Duval et al. 
(2009) for resolving the KDA for GRS molecular clouds. 

We utilize the procedure outlined in Simon et al. (2001) to convert from GRS 13 CO 
emission line parameters into an estimate of column density. This procedure assumes an 
excitation temperature of 10 K. As shown in Anderson et al. (2009), this excitation temper- 
ature is valid for the 13 CO emission surrounding young H II regions, and so is should also be 
valid for IRDCs. The conversion employs the equation 

iV(H 2 ) = 5.21 x 10 20 TAK, (1) 

where T is the line intensity and AV is the line width. The column density found using this 
procedure is accurate to within a factor of ~ 2 due to uncertainties in the conversion from 
13 CO to H 2 . 

To find the molecular column densities for both on- and off-cloud directions, we calcu- 
late the column density of all significant 13 CO molecular clumps in the X MM- Newton field of 
view by first fitting Gaussians to the average on- and off- cloud spectra. We transform the 
line parameters derived from these Gaussian into column densities using Equation 1. The 
application of Equation 1 computes the total column density along each line of sight, but 
give no information about whether the components lie at the near or the far distance. By 
examining the average 13 CO and H I spectrum of each molecular dump in the XMM-Newton 
field of view for H I SA, we determine whether the near or the far distance is appropriate 
for each molecular clump. We assign the column densities of each molecular clump calcu- 
lated from the average on- and off-cloud 13 CO spectra accordingly to the foreground and 
background column density numbers. 
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3.2. Atomic Column Density 

We use a similar procedure to calculate the Hi column density using the VGPS. The 
conversion to H I column density employs the equation 

iV(HI) = 1.82 x 10 18 T s J t(v) dv, (2) 

where T s is the Hi spin temperature, and v is given in kms -1 . For Hi, the spin temperature 
is much less well constrained compared to the excitation temperature of 13 CO. A good 
average value between the hot and cool Hi components is 150 K. Equation 2, requires first 
an estimate of the optical depth for the H I line using the equation 

THI = ~ ln ( l ‘ 5T=W) ' (3) 

where r H i is the optical depth of H I, X L is the H I line intensity, T s is the H I spin temperature, 
and 7 rg is the intensity from background sources. For simplicity, we assume no background 
emission. Assuming a spin temperature of 40 K lowers our calculated column densities by a 
factor of 2. 

For H I there is little information about whether the gas originates in front of or behind 
the cloud. As we did for 13 CO, we compute average on - and off-cloud spectra. A typical 
Hi spectrum does not show distinct, clean lines. We therefore compute the area under 
the curve directly, instead of fitting Gaussians, and convert this to H I column density using 
Equations 2 and 3. We split the column density distribution into foreground and background 
components, weighted by the distance to the cloud. Our values are in very good agreement 
with those in the the Leiden/Argentine/Bonn Hi survey (Kalberla et al. 2005) and the Hi 
survey of Dickey & Lockman (1990) as found using the HEASARC online tool. 3 

Figure 3 shows the average on- (solid curves) and off-cloud (dashed curves) spectrum 
of 13 CO (black curves) and Hi (gray curves) for Figure 3. For G2S.37+0.07, all 13 CO 
components with velocities less than the cloud velocity of 78.6 kms -1 lie at the near distance 
and correspond to absorption in the H I spectrum (this is much clearer when analyzing each 
clump individually). Also evident in this figure is the strong Hi absorption at the cloud 
velocity. The H I emission is less strongly peaked than the 13 CO emission. 

The results of our column density modelling are summarized in Table 2. This table 
gives the column density of atomic hydrogen, N( H), molecular hydrogen JV(H 2 ), and the 


3 http: / /hcasarc.gsfc.nasa.gov /cgi-bin/Tools/ w3nh/ w3nh.pl 



total column density, At = A r (H) + 2A(H2) in units of 10 22 cm for G28.3 i I 0.07 and 
G36.67-0.ll. The contributions to the column density are given for the foreground and 
background components ( “F” and “B” ) for the on— and off— cloud directions. The background 
component includes the contributions from the clouds themselves. 


4. Spectral Analysis 

Using our cleaned events files, we extract the X-ray spectra for the both on- and off- 
cloud regions in preparation for spectral fitting. Our on- and off-cloud regions are defined as 
before and correspond to high infrared extinction and high 13 CO column density. We extract 
spectra from both MOS detectors and create auxiliary response files (ARFs), redistribution 
matrix files (RMFs), and background spectra. We increase the signal to noise by binning 
these spectra so there are at least 20 counts per spectral channel. Below we describe the 
various model components we use to model the X-ray background and to separate foreground 
from background emission. We use XSPEC 4 Vll.3.2 (Arnaud et al. 1996) for the spectral 
analysis. 

There are two main sources of contamination in our spectra. First, there is residual soft 
proton contamination not removed through time filtering. This component should have a 
power law distribution of energies, so we model its emission using the XSPEC model pow/b. 
Secondly, there are two fluorescent instrumental lines of Al K a and Si Ka near 1.6 keV. We 
assume a Gaussian shape for these lines and fit them explicitly. 

Our models include three components for the X-Ray background, as in Snowden et al. 
(2000). Along each line of sight, first there is ~ 0.1 keV emission from the LHB. We fit 
the the LHB emission with the XSPEC model apec (Smith & Cox 2001), which represents 
collisionally-ionized diffuse plasma. Next, we assume that there are two components of the 
X-Ray background emission in the Galactic plane: a soft (< 1 keV) and a hard (few keV) 
component. For these components we use the XSPEC model wabs x apec, which represents 
collisionallv-ionized diffuse plasma absorbed by the Galactic column. The soft component 
itself can be divided into two sources: emission in front of the IRDC and emission from 
beyond the IRDC. Finally, there is extragalactic emission from unresolved active galactic 
nuclei (AGN) in the field. We model the AGN component using a powerlaw model with 
absorption, wabs x pow. where the powerlaw index is set to a = 1.46 (Chen et al. 1997). 
This is expected to be the dominant component above 1 keV (Lurnb et al. 2002). 


4 See http: / /xspec.gsfc.nasa.gov/docs/xanadu/ xspec 
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We show the model components and the column densities affecting them graphically in 
Figure 4 for G28.37+0.07. In this figure, molecular gas is displayed in black and atomic gas 
is displayed in gray. The white wedge is an exaggerated representation of the XMM-Newton 
field of view. Also shown on the figure are the location of IRDC G28.37+0.07, the local hot 
bubble (LHB), and the direction to the Galactic center (GC). The LHB size in the figure is 
exaggerated for clarity. Figure 4 shows that the AGN X-ray emission is affected by the total 
Galactic column. 

During the spectral fitting, we fit the spectra for each cloud simultaneously to improve 
the fit quality. We fix the abundance to Solar and use the Solar relative metal abundance of 
Anders & Grevesse (1989). We have five spectra for each cloud: on- and off-cloud spectra 
for both MOS detectors, as well as a ROSAT All Sky Survey (RASS Snowden et al. 1997) 
spectrum derived from the HEASARC X-ray Background Tool. 5 The RASS spectrum helps 
constrain the cosmic background at lower energies where the XMM-Newton detectors are less 
sensitive. We extract the RASS spectra from regions centered on our clouds of radius 0.5° and 
normalize the spectra to 1 arcminute 2 . We link the normalizations and temperatures of each 
spectral component for all five spectra to constrain the fits. We also link the temperature of 
the soft foreground component to that of the soft background component. 

The results of our spectral fitting procedure are shown in Figure 5 for the MOS1 detector. 
In Figure 5, the top curve is the on-cloud spectrum, and the lower curve is the off-cloud 
spectrum. The bottom curve is the RASS spectrum. The two prominent emission lines are 
the fluorescent instrumental lines of Al Ka and Si Ka. The left panel of Figure 5 shows the 
spectra and model for G28.37+0.07, and the right panel shows the same for G36.67— 0.11. 
It is evident from this figure that our model adequately fits the data. Also apparent from 
Figure 5 is the decreased counts for G36.67—0.11 due to the soft proton contamination. This 
results in greater uncertainty in the spectral fits. 

For G28.37+0.07, we run our fits in two trials, once with the column densities fixed to 
their calculated values (see §3) and once with the column densities allowed to vary (hereafter 
the “fixed” and “free” trials). We allow the column densities to vary up to 100% from 
our calculated values to account for uncertainties in converting line parameters to column 
densities. Due to insufficient counts, we are unable to allow the column densities to vary for 
G36.67— 0.11, and thus we run the spectral fits only with fixed values. 


5 ht t p : / / r osat . gsfc . nasa, gov / cgi-bin /Tools / xray bg/xray bg .pi 
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4.1. Spectral Fit Results 

The results of our spectral fits are shown in Figure 6 and listed in Table 3. Figure 6 
shows the best fit model components for G28.37+0.07 (left panels) and G36.67— 0.11 (right 
panels). The top panels show the models for the on-cloud directions, and the bottom panels 
show the same for the off-cloud directions. Table 3 lists the fitted temperature and emission 
measure (EM = J n 2 dl) for our four spectral model components, the reduced x 2 value, and 
the number of degrees of freedom in the fit. Errors in the table represent 90% confidence 
levels. Below we discuss the model components individually. 


4.1.1. Column Densities 

The column densities for G28.37+0.07 fit to similar, but generally lower, values com- 
pared to our calculated values (see §3). For the on— cloud spectra, the column density fits to 
1 2+°-03 x ^ q 22 cm -2 f or j} 1( , foreground component (we estimated 1.3 x 1 0 22 cm 2 ). The nor- 
malization for background column density fits to zero when the column densities are allowed 
to vary (see §4.1.4 below). Thus, the column density for this component is unconstrained by 
the spectral fits. The total column density affecting the Galactic ridge and AGN emission 
fits to 3.9ios x 10 22 cm" 2 for the on-cloud spectra (we estimated 5.2 x 10 22 cm" 2 ). For the 
off-cloud spectra, the column density fits to l.llgmt * cm ' 2 foreground to the cloud 
and to 2.91 q'2 x 10 22 cm- 2 for the total column density (we estimated 0.9 x 10 22 cm 2 and 
4.2 x 10 22 cm" 2 , respectively. 

Perhaps the most important column density parameter - the difference in column den- 
sity between the on- and off-cloud regions - is the same in the fixed and free trials for 
G28.37+0.07: 1.0 x 10 22 cm" 2 . The other discrepancies are discouraging, although the dif- 
ferences have little affect on the other model parameters. The largest effect can be seen in 
the emission measure of the foreground soft component (see below). Despite our care taken 
in estimating the column densities, errors in the conversion from 13 CO line parameters to 
proton column density are large and currently unavoidable. Longer or more sensitive obser- 
vations are required to more accurately disentangle the competing effects of absorption from 
the Galactic column and emission from the X-ray background. 


4.1.2. LHB Component 

The fixed and free column density trials for G28.37+0.07 both result in a LHB compo- 
nent temperature kT = 0.11 keV (logT/K - 6.11). For G36.67-0.ll. the LHB component 
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fits to a temperature of 0.089 keV(logT/K = 6.02. These temperatures for the LHB compo- 
nent are in good agreement with that found by many previous studies (see Snowden et al. 
1998, 2000; Kuntz k Snowden 2000; Henley, Shelton, k Kuntz 2007). The emission measure 
for the LHB component is less than what has been found by previous authors. We calculate 
an emission measure of 0.00055 cm -6 pc for the fixed trial and 0.00054 cm -6 pc for the free 
trial for G28.37+0.07. For G36.67-0.ll, it is slightly higher: 0.00063 cm- 6 pc. These values 
are three times less than that found by Henley, Shelton, k Kuntz (2007), ~ 10 times less 
than that found by Smith et al. (2007). This emission measure is in large part determined 
from the RASS 1/4 keV band data, and as the directions lie in the Galactic plane the surface 
brightness, and therefore the emission measure of the RASS data are relatively low. 


4-1.3. Soft Foreground Component 

The soft component fits to a temperature of 0.29 keV (logT /K = 6.53) in the fixed trial 
and to 0.30 keV (logT/K - 6.54) in the free trial for G28.37+0.07. For G36.67-0.ll, the 
soft component is not very well constrained. In fact, the only constraint we can place on 
this component is that it is less than 0.35 keV(logT/K = 6.61). Kuntz k Snowden (2008) 
find a similar temperature for this component: 0.24 keV. 


4 .I. 4 . Soft Background Component 

For the fixed trial of G28.37+0.07, the soft background component has little effect on 
the model. In fact, its removal has no impact on the quality of the fits, nor any impact 
on the value of the other parameters. The large column density in front of the background 
component causes significant absorption. The negligible effect this component has on our 
total model can be seen in Figure 6. For the free trial, the normalization for this compo- 
nent fits to zero. For G36.67-0.ll, the soft foreground component and soft background 
component have similar intensities. The difference between the fits seen in Figure 6 for the 
two clouds is due partly to the decreased column density to G 36.67— 0.11. Perhaps a larger 
factor, however, could be the greater uncertainty in the spectral fits for G36.67— 0.11 due to 
soft proton contamination. 
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4-1.5. Galactic Ridge Component 

The temperature for the warm Galactic ridge component fits to 1.61 keV (logT/K = 
7.27) in the free trial and to 2.34 keV (logT/K = 7.43) in the fixed trial for G28.37+0.07. 
This is a significantly higher than what has been found by previous authors. Kuntz k 
Snowden (2008) find a temperature of 0.71 keV for this component, with rather large error 
bars. Because of the decreased counts in the observation of G36.67-0.ll, we were forced to 
fix the Galactic ridge component to the temperature found for G28.37+0.07. When allowed 
to vary, this component would fit to an extremely high temperature (> lOkeV), which we 
deem unphysical, especially considering the results for G28.37+0.07. 

Figure 6 shows that for both clouds the Galactic ridge emission is the strongest compo- 
nent of the model between 2 and 5 keV. To determine if this component is required by the 
spectral fits, we remove the hard component from the models for G28.37+0.07. Its exclu- 
sion results in significantly worse fits: the reduced y 2 value increases from 1.11 to 1.41. We 
conclude that this is a necessary component of our spectral fits although its temperature is 
not very well constrained. 


4.2. O VII and O VIII Emission 

Oxygen is the most abundant line-emitting element at the temperatures of the LHB. In 
a million degree plasma, oxygen is primarily in the 0+ 6 state (O VIl), and the strongest lines 
detectable by XMM-Newton are near 0.5 keV. The emission from other charge states (most 
important for our purposes is O VIIl) can be used as a diagnostic to constrain the temperature 
of the LHB and verify the temperature found previously through spectral fitting. There are 
three lines of O vn that are all near 0.57 keV, and there is a Lya transition of O VIII at 0.65 
keV, 

We measure the intensities of the Ovil and Ovm emission from the local bubble by 
replacing the apec model component with a vapec model. The vapec model is identical to 
the apec model except that the atomic abundances may be individually modified. We set the 
oxygen abundance to zero and add two Gaussians representing the O vn and Ovm emission 
lines at 0.57 keV and 0.65 keV. For G28.37+0.07, when we fix all the model parameters to 
their previously fit values (except for the vapec temperature and normalization), we find 
that the Ovil line has an intensity of 3.41 j;q photons cm 2 s 1 sr 1 (hereafter line units, 
LU). The vapec temperature fits to a value of 0.098 keV. When we fix the vapec temperature 
to the LHB temperature found previously, 0.11 keV, we find an Ovil intensity of 3.3l 08 LU. 
We do not detect the Ovm line in either of these trials. The same procedure applied to 
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G36.67— 0.11 yields similar results. With the vapec temperature allowed to vary, the Ovil 
intensity is 1.9ti;| LU. The vapec temperature fits to 0.085 keV. When fixed to the 0.089 
keV value found previously, the Ovit intensity is 2.lt]J LU. We do not detect the Ovin 
line. 

Our values for the Ovil and Ovill lines are in rough agreement with what has been 
found by other authors. In a shadowing experiment to a nearby (~ 200 pc) filament, Hen- 
ley, Shelton, & Kuntz (2007) found 3 LU for the Ovil line and did not detect the 
Ovill line. Kuntz & Snowden (2008) found 1.75+/-0.7 LU for the Ovil emission at (l,b) 
— (lll.°14, l.°ll). In a shadowing experiment of MBM 12, a high latitude cloud ~ 50 to 
~ 300 pc distant, Smith et al. (2005) found 1.79 ± 0.55 LU for O vn and 2.34 ± 0.36 LU for 
O VIII. They note that this unusually strong Ovill emission can be explained by solar wind 
charge exchange (SWCX) contamination in their spectra. The SWCX emission is caused 
by electron transitions between neutral atoms and highly ionized species in the solar wind. 
It has a potentially large effect in X-ray data, especially following a coronal mass ejection 
event. Snowden, Collier & Kuntz (2004) showed the large effect on the oxygen lines that 
SWCX emission can have. In a followup study of MBM12 with Suzaku , Smith et al. (2007) 
found 2.34 ± 0.33 LU for Ovil and 0.77 ± 0.16 LU for Ovill line emission. Using the X-ray 
Quantum Calorimeter aboard a sounding rocket, McCammon et al. (2002) found 4.8 ± 0.8 
LU for the Ovil emission and 1.6 ± 0.4 LU for the O VIII emission. 

A consensus appears to be emerging that most of the Ovil emission arises from the 
LHB (with a contribution from SWCX). Strong O VIII emission, however, cannot arise from 
the LHB as this would imply a temperature inconsistent with other determinations. This 
emission must arise from SWCX (Snowden, Collier &: Kuntz 2004) or from a cool compo- 
nent distributed throughout the Galaxy (Kuntz & Snowden 2008), likely with contributions 
from both sources. Our results are in agreement with this interpretation, although the cool 
component in our models does not produce appreciable intensity in the O VII line - our data 
are consistent with all the Ovil line intensity originating from the LHB. That we do not 
detect the O VIII line is not surprising given the large absorbing column and low expected 
line strength. 


5. Image Analysis 

We produce exposure corrected, background subtracted X-ray images to verify the re- 
sults from the spectral analysis, and to search for a morphological match between the molec- 
ular emission and an X-ray decrement. Using the XMM-ESAS software, we create these 
images in three energy bands: 0.35-1.25 keV, 1.25-2.00 keV, and 2-10 keV. We model two 



- 14 


sources of background in our images: the particle background and residual soft proton con- 
tamination. These background components must be removed as cleanly and completely as 
possible to ensure accurate results. 

Using the filter-wheel-closed data, the XMM-ESAS software models the particle back- 
ground (see Kuntz & Snowden 2008). These data are dominated by the instrumental back- 
ground because the chips were not exposed to the sky. We do not include blank sky data 
in the modelling of either the particle background nor the soft proton contamination. The 
blank sky files may suffer from residual soft proton contamination and solar wind charge 
exchange. Additionally, the blank sky files contain the cosmic background which we are 
attempting to observe. Using these data would not only add error to our models, but would 
also remove the signal of interest. 

The XMM-ESAS software also allows for the creation of a model image of any residual 
soft proton contamination. Kuntz & Snowden (2008) characterize the reasonable ranges 
of the variations in both the spectrum and in the spatial distribution over the detectors. 
XMM-ESAS incorporates this characterization in the creation of a model soft proton image. 
The energy distribution of the soft proton contamination is determined from the power law 
index and normalization found in the spectral fits. For this input, we use the values found 
in the spectral fitting. 

We subtract the model background and soft proton images, combine the data from 
the two MOS detectors and smooth the resultant image. The result of these operations (an 
exposure corrected, background subtracted image) is shown in Figure 7 for the 0.35-1.25 keV 
band. In Figure 7, the left panel is of G28.37+0.07 and the right panel is of G36.67-0.ll. 
The strongest strongest X-ray shadow should be present in the 0.35-1.25 keV image. In the 
images of both clouds there is some indication of a shadow, but it not the clearly defined 
shadow one would expect, nor does it have the same depth one would expect. Furthermore, 
the same shadow that appears in the 0.35-1.25 band also appears in the 2. 0-8.0 keV band. 
As there should not be strong absorption in the 2.0-8. 0 keV band, this indicates that the 
“shadow” is a product of the detector, and not of the cloud itself. We conclude that the 
imaging analysis supports the spectral analysis and finds no X-ray shadow of the cooler 
component to the X-ray background. 


6. Summary 

We have conducted X-ray shadowing experiments on two infrared dark clouds near 
30° at distances of 2 and 5 kpc from the Sun to determine the spatial distribution 
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of the X-ray background. We used Hi and 13 C0 data to calculate the total atomic and 
molecular column densities for both the on- and off-cloud directions. Both clouds have very 
high proton column densities of > 4 x 10 22 cm” 2 . This high column density should absorb 
nearly all the soft background X-ray flux. We find that the cool, diffuse X-ray background 
must originate forground to the clouds, within a few kpc of the Sun. Our results are further 
evidence that X-ray emitting plasma is distributed throughout the disk of our Galaxy (see 
Park et al. 1997; Kuntz & Snowden 2008). 

We found that the X-ray background is best fit with a three component model with 
contributions from the local hot bubble (LHB), a soft component, and a hot component 
from the Galactic ridge. The LHB emission is best fit with a value of ~ 0.1 keV (log T / K = 
6.06) for both clouds, in agreement with what has been found in previous studies. The soft 
component is best fit with a temperature of ~ 0.3 keV (log T /K = 6.54). This is the dominant 
source of emission between 0.7 and 1.0 keV. The Galactic ridge component is best fit at values 
near 2 keV. This is the dominant source of emission in our spectral models between 2 and 
5 keV. Our spectral fits show that all significant emission below 1 keV can be attributed to 
the LHB. Image analysis including accurate modelling of any background components failed 
to reveal an X-ray shadow, supporting this results of the spectral modelling. 

We detect the OviI line for both clouds at levels of ^ 3 photons cm” 2 s” 1 sr” 1 . The 
O VII line intensity is consistent with this emission originating from the LHB. These values 
are roughly in agreement with what has been found in previous studies. We do not detect 
the O VIII line for either cloud. 


This publication was based on observations obtained with XMM-Newton , an ESA sci- 
ence mission with instruments and contributions directly funded by ESA Member States and 
NASA. This work was supported by NASA XMM- Newton Guest Observer grants including 
NNX06AG73G and NNG05GP68G. This publication makes use of molecular line data from 
the Boston University-FCRAO Galactic Ring Survey (GRS). The GRS is a joint project of 
Boston University and the Five College Radio Astronomy Observatory (FCRAO), funded by 
the National Science Foundation under grants AST-9800334, AST-0098562, AST-0100793, 
AST-0228993, & AST-0507657. This publication makes use of data from the VLA Galactic 
plane survey (VGPS). The National Radio Astronomy Observatory is a facility of the Na- 
tional Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 
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Fig. 1. — IRDCs shown in absorption against GLIMPSE Spin emission. The left panel shows 
G28. 37+0.07 and the right panel shows G36.67-0.ll. The contours are 13 CO integrated 
intensity from the GRS. Tick-marks on the contours point downhill, toward decreasing 13 CO 
integrated intensity. The large black circles represent the XMM-Newton field of view (~ 30'). 
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Fig. 4. — Components of the X-ray background attenuated by the Galactic column for 
G28.37+0.07. The column density Molecular gas from the GRS is shown black while atomic 
gas from the VGPS is shown in gray. The Local Hot Bubble (labeled LHB), the direction of 
our observation (white wedge) and the direction towards the Galactic center are also shown. 
The size of the wedge and the LHB are slightly exaggerated for clarity. IRDC G28.37+0.07 
is associated with a large dense molecular cloud 5 kpc from the Sun.) 
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Fig. 5. — Spectral data and models for G28.37+0.07 (left panel) and G36.67— 0.11 (right 
panel) from the MOS1 detector. Shown are the on-cloud (black), off-cloud (dark gray), and 
RASS (light gray) spectra. The bottom section of each panel shows the fit residuals. 
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Fjg g. — Model spectral components from the XSPEC fit to the on-cloud (top panels) 
and off— cloud (bottom panels) spectra for G28.37+0.07 (left. panels) and G36.6/—0.11 (right 
panels) from the MOS1 detector. 
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Fig. 7. — Smoothed XMM-Newton images in the energy range 0.35 keV to 1.25 keV of 
G28. 37+0.07 (left panel) and G36.67-0.11(right panel). The contours are the same GRS 
contours shown in Figure 1. There is no obvious X-ray decrement at the location of the 
molecular clouds. 




Table 1. Cloud Properties 


Name 

I 

b 

D 

Size 

N(H 2 ) 


° 

o 

kpc 

sq. ' 

cm” 2 

G28.37+0.07 

28.37 

+0.07 

5.0 

42 

21 

G36.67~~0.il 

36.67 

-0.11 

2.0 

19 

19 


Table 2. Cloud Column Densities^ 



G28.3' 

ON 

r+0.07 

OFF 

G36.67- 

ON 

~0.11 

OFF 

F 

B 

F 

B 

F 

B 

F 

B 

W(H) 

0.5 

1.0 

0.5 

1.0 

0.5 

1.0 

0.5 

1.0 

iV(H 2 ) 

0.4 

2.1 

0.2 

1.6 

0.1 

1.4 

0.1 

0.9 

N r 

1.3 

5.2 

0.9 

4.2 

0.7 

3.8 

0.7 

2.8 


+ A1J column densities are in units of 10 22 cm' 2 . 
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Table 3. Spectral Model Parameters 




G28.37+0.07 


G36.67- 

-0.11 


Column Fixed 

Column Free 

Column Fixed 

Parameter 

Value 

Error 

Value 

Error 

Value 

Error 

LHB 

kT(keV) 

EM(cm“ 6 pc) 

0.11 

0.00055 

+0.0037 
-0.0093 
+7.6E — 5 
- 7.0E — 5 

0.11 

0.00054 

+0.00830 
-0.00097 
+6E — 5 
— 7.1E-5 

D.089 

0.00063 

+0.0047 

-0.0092 

+3.6E-5 

-6.8E-5 

Cool Foreground 
kT(keV) 

EM (cm 6 pc) 

0.29 

0.046 

+0-0150 

0.0076 

+0.0069 

0.0098 

0.30 

0.034 

+0.028 

0.015 

+0.0025 

0.0100 

< 0.35 
0.0023 

+0.00037 

0.00150 

Cool Background 
kT(keV) 

EM (cm “ 6 pc) 

0.29 

0.0019 

+0.157 

-0.0019 

0.30 

0.00 

+3.9E-5 

- 0.000 

< 0.35 
0.058 

+0.0190 

-0.0087 

Galactic Ridge 
kT(keV) 

EM (cm ' 6 pc) 

x 2 

dof 

1.61 

0.022 

1373.5 

1236 

+0.25 

-0.14 

+0.0015 

0.0047 

2.34 

0.012 

1348.1 

1231 

+ 0.35 
-0.28 
+0.0021 
0.0011 

1.61 

0.0082 

352.4 

300 

+0.0016 

0.0026 



